! ***************************************
finish
/clear
/filn, e341
/title, mechanical stress analysis of skirt supporting zone of hydrogenation reactor
/units,si
! ********* 趨*********
pi=3.1415926                            ! Բ
tongt_di=2813                           ! Ͳھ
tongt_t=87                              ! Ͳں
tongq_h=615                             ! Ͳ嵽ȹ֧дĸ߶
fengt_t=52                              ! ͷں
fengt_di=2833                           ! ͷھ 
xyhr1=20                                ! ڲԲ뾶
xyhr2=20                                ! Բ뾶
h1=1500                                 ! ȹ߶
h2=1000                                 ! Ͳ峤
xqunz_h=70                              ! hζͼȹӲֱ߸߶
qunz_t=22                               ! ȹں                                 
fengt_do=fengt_di+2*fengt_t             ! ͷ⾶
presi=8.83                              ! ѹ
presa=presi*(tongt_di**2)/((tongt_di+2*tongt_t)**2-tongt_di**2)  
	                                ! ѹͲ˲ƽ
mass=270000                             ! 豸
presg=4*mass*9.8/(pi*((tongt_di+2*tongt_t)**2-tongt_di**2)) 
                                        ! 豸Ͳ˲Ӧ 
EXX1=2e5                                ! ϵĵģ
NUXY1=0.3                               ! ϵĲɱ
!****************ǰ***************************
/prep7
et,1,plane82                            ! 嵥Ԫ
mp,ex,1, EXX1                           ! ϵĵģ
mp,nuxy,1, NUXY1                        ! ϵĲɱ
keyopt,1,3,1                            ! 趨Գѡ
!****************ģ***************************
k,1,0,0 	                        ! ĵ
k,2,tongt_di/2+tongt_t,0                ! hζͼ϶
k,3,tongt_di/2,0                        ! hζͼ϶ڲ
k,5,tongt_di/2+tongt_t,-tongq_h         ! hζͼȹӲ
k,6,tongt_di/2+tongt_t-qunz_t,-tongq_h	! hζͼȹӲڲ
k,7,kx(6)-xyhr1,-(tongq_h-xqunz_h)	! hζͼȹӲԲ
l,2,3                                   ! 
l,2,5
l,5,6

circle,7,xyhr1                          ! ȹӲԲԲ
ldele,5,7,1,1                           ! ɾԲ

l,6,4                                   ! 
k,7,0,ky(8)                             ! ¶ؼ7ΪԲǸߵĵ
l,7,8                                   ! hζͼԲߵĸ
circle,1,fengt_do/2                     ! ڶӦԲ
ldele,7,9,1,1                           ! ɾԲ

lsbl,10,6                               ! ߶
ldele,7,,,1                             ! ɾ߶

l,8,10                                  ! ٴԲߵĸ
k,30,kx(10)-xyhr2,ky(8)-xyhr2           ! 帨ؼ
k,40,kx(8),ky(30)                       ! 帨ؼ
 	
l,40,30                                 ! ɸȷڲԲ
circle,1,fengt_do/2+xyhr2               ! ԲȷڲԲ


lsbl,7,12                               ! 󽻵ȷڲԲ
ldele,13,14                             ! ɾ߶
kdele,30                                ! ɾĹؼ
ldele,9,11,1,1                          ! ɾ߶

circle,14,xyhr2                         ! ڲԲԲ
ldele,7,,,1                             ! ɾĹؼ߶
ldele,10,11,,1                          ! ɾĹؼ߶

l,1,14                                  ! ɸ
l,9,14                                  ! ɸ
k,5000,kx(14),ky(14)                    ! ٶڲԲ
lsbl,9,7                                ! ߶
ldele,12,,,1                            ! ɾ߶

l,7,14                                  ! ɸ
lsbl,6,10                               ! ߶
ldele,12,,,1                            ! ɾ߶
lsbl,8,7                                ! ߶
ldele,6,,,1                             ! ɾ߶

k,1,0,0                                 ! ¶
circle,1,fengt_di/2                     ! ڱڵԲ
ldele,6,8,1,1                           ! ɾ߶

k,40,kx(3),ky(7)                        ! 帨ؼ
l,3,40                                  ! ɸ
l,12,14                                 ! ɸ
lsbl,12,6                               ! ߶
ldele,8,,,1                             ! ɾ߶

ldele,3                                 ! ɾ߶
k,100,kx(5),ky(5)-h1                    ! ȹ߶εĹؼ
k,200,kx(6),ky(6)-h1                    ! ȹ߶εĹؼ
l,100,200                               ! ȹ߶
k,300,kx(3),ky(3)+h2                    ! Ͳ߶εĹؼ
k,400,kx(2),ky(2)+h2                    ! Ͳ߶εĹؼ
ldele,1                                 ! ɾ߶
l,300,400                               ! Ͳ߶ 
l,300,11                                ! Ͳڱ߶
l,400,2                                 ! Ͳ߶
l,5,100                                 ! ȹ߶
l,6,200                                 ! ȹڱ߶
k,800,kx(2),ky(11)                      ! Ӧ·ؼ
k,900,kx(5),ky(4)                       ! Ӧ·ؼ
k,6000,kx(300),ky(300)-400              ! Ӧ·ؼ
k,7000,kx(400),ky(6000)                 ! Ӧ·ؼ
l,1,5000                                ! ɸ
lsbl,13,15                              ! ߶
ldele,2                                 ! ɾ߶   
a,12,14,10,11,300,400,2,800,900,5,100,200,6,4,8,9,7     ! 

!***************************************
aesize,1,10                             ! 嵥Ԫߴ
mshkey,0                                ! ʷ
amesh,1                                 ! ʷ
lrefine,4,,,1                            ! Բ
lrefine,11,,,1                           ! ڲԲ
fini 
!*******************************************
/solu
sfl,6,pres,presi                        ! ڱʩѹ
sfl,16,pres,presi                       ! ڱʩѹ
sfl,17,pres,presi                       ! ڱʩѹ
sfl,1,pres,presg-presa                  ! Ͳ˲ʩƽغ
dl,3,1,uy,0                             ! ȹ׶߶ʩλԼ
dl,7,1,ux,0                             ! ͷԳʩxλԼ
solve                                   !          
finish



!*******************************************
/post1
plnsol,s,int,0,1                        ! Ӧǿͼ
!****************Ӧ·1--1****************
node11=node(kx(800),ky(800),kz(800))    ! ڽڵ
node12=node(kx(11),ky(11),kz(11))
path,set1,2                             ! ·Ƽ·       
ppath,1,node11                          ! ݽڵ㶨·
ppath,2,node12                          ! ݽڵ㶨·
prsect,-1                               ! г·Ի
!****************Ӧ·2--2****************
node21=node(kx(4),ky(4),kz(4))
node22=node(kx(900),ky(900),kz(900))
path,set2,2        
ppath,1,node21
ppath,2,node22
prsect,-1
!****************Ӧ·3--3****************
node31=node(kx(10),ky(10),kz(10))
node32=node(kx(7),ky(7),kz(7))
path,set3,2        
ppath,1,node31
ppath,2,node32  
prsect,fengt_di/2
fini
